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Abstract 

The IceCube Observatory at the South Pole is composed of a cubic kilometer scale neutrino telescope 
■ buried beneath the icecap and a square-kilometer surface water Cherenkov tank detector array known 

as IceTop. The combination of the surface array with the in-ice detector allows the dominantly elec- 
' tromagnetic signal of air showers at the surface and their high-energy muon signal in the ice to be 

04 \ measured in coincidence. This ratio is known to carry information about the nuclear composition 

of the primary cosmic rays. This paper reviews the recent results from cosmic-ray measurements 
performed with IceTop/IceCube: energy spectrum, mass composition, anisotropy, search for PeV 7 
sources, detection of high energy muons to probe the initial stages of the air shower development, and 
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ly-^ ■ study of transient events using IceTop in scaler mode. 
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■ I. Introduction 

For 100 years since their discovery, primary cosmic rays have been measured with different 
techniques and at different locations on Earth[lJ. Particularly interesting are those cosmic 
particles in the energy range between 3 -lO^"^ eV (300 TeV or 0.3 PeV) and lO^^ eV (1000 EeV), 
O \ whose origin, composition, and energy spectrum remain not fully understood. Cosmic mag- 

netic fields permeating all space prevent the localization of the sources that produced the 
^ , charged particles observed at Earth. In addition, the relatively low flux of these primaries 

(about 1 m^^ yr^^ at 10^^ eV) does not allow for direct measurements. Subsequent interac- 
tions in the atmosphere result in air showers of secondary particles that are sampled at groimd. 
^ ■ These secondaries also release Cherenkov and fluorescence light detected with telescopes. 

. The IceCube Observatory is a three-dimensional cosmic-ray air shower detector (Fig. [1). 

The surface component, IceTop, is an array of water Cherenkov tanks that samples the shower 
at the surface. The deep detectors of IceCube (between 1.45 km and 2.45 km below the surface) 
measure the signal from penetrating muons, which have 500 GeV or greater at production 
in the atmosphere. Events seen in coincidence by both the surface and the deep detectors 
CS| ■ (coincident events) offer a unique view of cosmic-ray air showers because they carry informa- 

tion from two different regions of shower development. The surface shower is dominated by 
low-energy photons, electrons, and muons produced throughout the cascade, while the high- 
energy muons reflect the early stage of shower development. In addition to coincident events, 
it is also possible to use IceTop as a stand-alone air shower array and to use the deep detector 
d \ by itself as a muon detector. Together, the three types of events allow measurement of the 

energy spectrum, composition, and anisotropy of the primary cosmic rays from about 10 TeV 
to about 1 EeV. 

The main goal of IceCube is to detect and measure high-energy neutrinos of extra-terrestrial 
origin, as described in the companion review article in this journall2|. The high-energy neutri- 
nos are expected to point back to the cosmic-ray sources in which they are produced|3|. From 
the point of view of neutrino astronomy, IceTop serves as a partial veto against atmospheric 
background in the deep detector. 

The emphasis of this paper is to review the apparatus and its performances, and the recent 
physics results of IceCube as a cosmic-ray detector. The performances will be reviewed in 
Sec. 121 Current results from cosmic-ray mass composition and energy spectrum analysis will 
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Figure 1: Sketch of the IceCube Observatory labeled with some of its main features. The different colors at the 
surface identify different deployment stages of the detector. IceCube in its 2006-07 configuration is 
shown in red and termed as IT26/IC22 (26 IceTop stations/22 in-ice strings). Other configurations 
are IT40/IC40 (2007-08) in green, 1T59/IC59 (2008-09) in violet, 1T73/IC79 (2009-10) in blue, and 
IT81/IC86 (2010-11) in yellow. 



be discussed in Sec.|3] The anisotropy measured with IceTop data is presented in SeclH Finally, 
this paper will also report on recent results including PeV j search (Sec. |5j, high momentum 
in-ice laterally separated muons (Sec. [6]l, and detection of transient events (solar flares and 
gamma-ray bursts) with IceTop (Sec. 01. 

11. Cosmic-ray Detection and Reconstruction 

The IceTop detector[4] is a surface array of 162 cylindrical Cherenkov tanks installed at the 
2835 m altitude of the South Pole surface (atmospheric depth of about 680 g/cm^). The tanks 
are 1.86 m in diameter and 1.10 m in height. Filled with clear ice to a depth of 0.90 m, they 
operate on the same principle as the water tanks of the Haverah Park experiment|51 and the 
Pierre Auger Observatory[6J. To minimize accumulation of drifting snow, the tanks have their 
top surface level with the surroimding snow. Nevertheless, a variable and not negligible snow 
coverage is measured and accounted for when analyzing data. 

Pairs of tanks, 10 m apart from each other, localize 81 stations that are distributed over an 
area of about 1 km'^ on a triangular grid with mean spacing of approximately 125 m. An in-fiU 
array for denser shower sampling is configured in the center of IceTop by using 3 stations at 
smaller distances together with their 5 neighboring stations. Each IceTop tank contains two 
standard IceCube digital optical modules] 7, 8| (DOMs), which include light sensor and read- 
out electronics. To enhance the d5mamic range, the DOMs of each tank are run at 2 different 
gains (low and high) with resulting effective thresholds of 20 and 200 photoelectrons (pe), re- 
spectively. Cosmic-ray muons, hitting the tanks at an approximate rate of 2 kHz, provide the 
basic calibration of tanks. The signal spectrum of a tank consists of a low-energy electromag- 
netic component and a muon peak at higher energy. The tank signals are calibrated in units 
of "vertical equivalent muon" or VEM whose definition is based on the muon peak of the tank 
spectrum. A vertically through-going muon of about 1 GeV produces approximately 125 pe in 
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Figure 2: (a) Coincident event recorded by IceTop/lceCube in 2010 (1T73/IC79). The triggered DOMs are indicated 
with colored spheres with radii proportional to the signal, (b) Lateral distribution of tank signals in VEM 
for the event in Fig. \2(a)\f itted to a double logarithmic parabola (see text for details), (c) Muon energy 
loss dE/dX as a function of the in-ice depth for the event in Fis. \2{a)\ The stochastic energy loss is about 
880 GeV/m above the fitted average energy loss (black line). The dashed line is the average energy loss 
after removal of the stochastic peaks. 

a high-gain DOM. In addition to measuring signal amplitudes, the IceTop DOMs can record 
the counting rate of low-energy cosmic rays (scaler mode). The rates are available for helio- 
spheric studies of solar modulation and transient events such as solar flares and gamma-ray 
bursts (GRBs). 

To initiate the readout of a DOM, its neighbor in the other tank at the same station has to 
be hit (hard local coincidence or HLC). This suppresses the background of accidental signals 
caused by isolated muons and allows a good angular resolution. The basic IceTop air shower 
trigger requires 6 DOMs to "launch" or report signals within a 5 ^s time window. Thus, even 
if only high-gain DOMs are above threshold, this trigger includes all 3-station events. The 
rate for these events is 30 Hz. Once the trigger is formed, additional signals from single tanks 
are added if they occur within a 20 f/s time window (soft local coincidence or SLC). Counting 
SLCs is useful to evaluate the shower muon content. These SLC signals are also read out when 
an rn-ice trigger forms either with or without the presence of HLC signals. In this case, SLCs 
can be used to veto dow-going cosmic rays when measuring up-going neutrinos. 

There are different topologies of events that are relevant for analysis of cosmic-ray data 
with IceTop /IceCube. This paper will concentrate on results obtained with analysis of events 
caused by nearly vertical (zenith angles 6 < 37° or cosO > 0.8) and contained air showers 
(Fig. |2(a)| . These showers are reconstructed with their axis crossing both parts of the detector. 
The effective area of IceCube for such coincident events is ApsO.15 km^sr. The maximum 
energy above which the intensity is too low to obtain enough events for analysis is about 1 EeV. 
The effective area rises to 0.4 km^sr for detecting events by IceTop alone. 
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The surface shower particle density decreases rapidly with the distance from the shower 
axis (lateral distribution function). This lateral distribution carries information about the en- 
ergy of the primary particles. The charge expectation value S in an IceTop tank at distance r 
from the shower axis is described with a "double logarithmic parabola" (Fig. |2(b)| as follows l9l 

where Ry^f = 125 m, Sy^f is the charge in VEM at Rygf, and k = 0.303. The parameter Sy^f 
is thus referred to as S125 and is a measurement of the shower size. The signals measured 
between about 30 m and 300 m from the shower axis are quite well described by Eq. [1] for 
primary energies in the range 1-100 PeV and arrival directions with zenith in the range 0°- 
40°. 

The arrival time behind the shower plane, Af (r), as a function of the lateral distance r from 
the shower axis is described by the sum of a parabola and a Gaussian function, both symmetric 
around the shower axis, with the following form 

"' + ^'fexpf^' 



Af(r) =flr^ + Mexp( ^) -1), (2) 



where a = 4.823-10~* ns/m^, h = 19 Al ns, and ay = 83.5 m. The energy, zenith angle, 
and mass dependence of a, h, and cXy is currently under investigation in simulations. When 
reconstructing the shower direction n, Eq.|2]accoimts for the shower front curvature as follows 

t{x) =tQ + ^^^n + M{r), (3) 

where t(x) is the tank signal time, x the tank position, Xc the shower core position, and the 
time when the core reaches the surface. For events triggering 5 or more stations, Eq. |3]is fitted 
to the measured signal times with 5 free parameters, 3 for core position and time, and 2 for 
shower direction. At the same time, Eq. [1] is fitted to the measured tank signals with 2 free 
parameters, S125 and j6. A maximum likelihood method is adopted to obtain the best fit to 
the measurements[4J. For events triggering 3 or 4 stations, only S125 and /5 are kept free after 
fitting the measured signal times to obtain core position and time, and shower direction. 

The position of IceTop at the high altitude of the South Pole makes it possible to sample 
ground particles near the shower maximum, thus reducing significantly the effects of fluctu- 
ations and allowing accurate measurements. Events of 5 or more stations are expected from 
primaries of energy between PeV and EeV, whereas 3 or 4 station events set the threshold for 
cosmic ray detection with IceTop to about 300 TeV. Lower energies (J^lOO TeV) are expected if 
the events triggering the denser in-fill array are considered. The energy resolution is estimated 
to be ^0.1 in units of log2Q(E/GeV) above about 1 PeV, reaching about 0.05 for E > 10 PeV. The 
angular resolution ranges from about 1° at 10^^ eV to about 0.2° at 10^'^ eV. The core resolution 
ranges from about 15.0 m at 10^^ eV to about 6 m at 10^"^ eV. 

Analysis of IceCube coincident events aims to give clear insights into the nuclear com- 
position of cosmic rays for energies that span from PeV to EeV. For a given primary energy, 
penetrating muons are more abundant in iron showers than in proton showers since their de- 
velopment starts higher in the atmosphere. The muon bundle size is therefore larger for iron 
showers. On the other hand, proton showers reach their maximum development deeper in the 
atmosphere than iron showers and the ratio of muons to electrons and photons is therefore 
smaller at the surface. The number of in-ice muons or muon multiplicity is closely related 
to the amount of energy deposited in the detector, which is proportional to the amount of 
Cherenkov light generated. For example, depending on the mass of the primary particle (pro- 
ton or iron), an event of 5-10^^ eV is expected to carry 30 to 80 muons with sufficient energy 



4 



Preprint • To Appear in Mod. Phys. Lett. A. 



to reach a depth of 1500 m and deposit 510 eV to 1510 eV in the deep detector. Analo- 
gous to surface measurements of tank signals, the in-ice DOM signals (in photoelectrons) can 
be described in terms of a lateral distribution, which is a function of the distance from the 
muon bundle track at a given slant depth from the surface. This function is dominated by a 
decaying exponential whose slope is the attenuation length of light in the ice|10|. The bundle 
size, termed as X70, is defined as the in-ice signal measured at a slant depth of 1950 m and at a 
perpendicular distance from the track of 70 m. This observable has been used in combination 
with S125 to discriminate between light and heavy primary masses (see Sec. |3]l. A different 
approach explored the reconstruction of the muon bimdle energy loss dE / dx to find compo- 
sition sensitive properties of bundles[llj. The muon bundle energy loss as a function of the 
slant depth is a convolution of the shower muon multiplicity, the muon energy distribution 
and the energy loss of a single muon. For iron showers, the bundle energy loss is expected to 
be greater. Furthermore, for the same amount of energy deposited, stochastic losses along the 
bundle track are larger in proton showers, which produce more high energy muons. 

The reconstruction of air showers with IceTop/IceCube is affected by several uncertainties. 
These uncertainties have been extensively investigated and several advances are expected in 
the future results. Pressure variations at the surface and snow coverage affect the measurement 
of S125 with the latter giving the largest contribution. For a given primary energy and arrival 
direction, higher pressure makes the primary interaction depth shallower thus reducing the 
shower size. The snow over each tank is physically measured every year and also estimated 
from the muon/electron ratio in calibration curves. The snow mainly affects the response 
of the tank to the electromagnetic component of the shower front thus affecting the energy 
threshold. The shower size is corrected for pressure and snow coverage. Other systematic 
imcertainties affect the in-ice reconstruction: the ice model used to describe the properties 
of the photon propagation^ (±10%), the DOM efficiency 1131 (±10%), and the muon rate 
seasonal variation ll4J . 

III. Cosmic-ray Energy Spectrum and Nuclear Composition 

Most galactic cosmic rays are believed to be accelerated in the blast waves (diffusive shocks) of 
nearby supernova remnants fTSlI (SNRs) and in some cases (extended sources /strong magnetic 
fields) can escape the acceleration region with energies up to about 10^^ eV. The signature of 
these sources is the gradual steepening of the cosmic-ray flux at a few 10^^ eV, called knee. 
The well-known power-law form of the cosmic-ray energy spectrum dN/dE oc E^^-'^ changes 
its spectral index to about -3.0. The knee is interpreted as the maximum energy reached 
through acceleration in SNRs. This energy scales with the charge of the nucleus. Therefore, 
the spectrum of galactic cosmic rays is predicted to end at energies of about 10^^ eV with 
the heaviest elements being accelerated] 16 1. Depending on whether the transition is from 
galactic iron to extragalactic proton or from galactic iron to extra-galactic mixed composition 
of different nuclei, the transition energy is predicted[17J at few 10^'' eV or few 10^^ eV. At 
energies between 10^^ eV to 10^'' eV, all air shower experiments observe energy dependent 
changes of composition that are compatible with an increase of the measured average mass 
of cosmic rays. Above 10^'' eV and up to 10^^ eV measurements of composition indicate a 
decrease of the cosmic-ray average mass| 18|. The cosmic-ray flux up to 10^'' eV is believed to be 
mostly dominated by the contribution of galactic sources. The "fingerprint" of the transition 
to the extra-galactic contribution is expected in the measurement of the mass composition 
at this energy and above. At the highest energies (above 10^^ eV), measurements from the 
Pierre Auger Observatory, and HiRes/Telescope Array give opposite results 1 19t l20ll2Tl . These 
experiments are extending their detectors to reach about 10^^ eV and below li22l l23l l24i . and 
their measurements will overlap with IceCube measurements. 

IceCube/IceTop allows for precise measurements of the primary energy spectrum in a wide 
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Figure 3: Energy spectra obtained with data oflceTop running in different configurations (1726, IT40, and 1T73). 
For the IT26 and IT40 spectrum, the systematic uncertainty is also shown. 



energy range that reaches the energy threshold of the largest air shower arrays and is sensitive 
to cosmic-ray nuclear composition changes. The first analysis to determine the all-particle 
energy spectrum with IceTop l|25l was based on data of IT26 (area of 0.094 km^) taken between 
June and October 2007. In this analysis, the measured shower size (S125) spectra in three zenith 
angle ranges up to 46° are unfolded or de-convoluted into the estimated energy spectrum. The 
spectrum obtained in the energy range between 1 and 100 PeV is shown in Fig. |3] along with 
other spectra measured with IceTop data and discussed in this section. The energy at which 
the lines with the minimum and maximum slope before and after the knee intersect identifies 
the knee position. Assuming pure iron, the energy spectra measured in different zenith angle 
ranges have been shown to disagree and prove that pure iron primaries can be excluded in 
the energy range up to 25 PeV. A consistent interpretation of the spectra measured at different 
zenith angles requires a mixed composition. For a two-component model, the knee measured 
in the IT26 spectrum is at 4.32 PeV and the spectral index above the knee is -3.11. An indication 
of a flattening of the spectrum above 22 PeV is also observed with a spectral index changing 
to -2.85. 

A preliminary measurement of the spectrum with IT73 has been obtained by analyzing 11 
months of data (June 15 to May 13, 2010)|26H2Zl|28l (Fig.|3ll. The statistics is nearly 40,000,000 
events between about 0.3 PeV and 1 EeV and for cos 9 > 0.8. Of these events, about 200 
are found above about 200 PeV. With enhanced precision, better energy resolution, and larger 
statistics above 100 PeV, this measurement confirms the earlier result of a flattening observed 
with IT26 data and reveals the spectrum structure in the steepening between the knee and ankle 
above about 100 PeV. It also appears that the spectrum is not well described by a single power 
law. All showers with 5 or more stations and with cos 6 > 0.8 are considered in this analysis. 
The energy spectrum, now measured between 1 PeV and 1 EeV, will be extended to lower 
energies, down to 300 TeV. This can be achieved by selecting small events in IceTop 1 29 1 and 
is particularly interesting in view of the recent measurements with ATIC ||30| and CREAM 1 3 1) 
balloon-borne calorimeters that are providing increased statistics and a new view of the region 
around 100 TeV/nucleus. 

In Fig. m the IceTop spectrum is compared to the spectra obtained with recent measure- 
ments of KASCADE-Grande (Karlsruhe, Germany, 110 m a.s.l.)||3l |35l, Tibet Array (at Tibet 
Yangbajing, 4300 m a.s.l.)|36. 37|, GAMMA (on the south side of Mount Aragats in Armenia, 
3200 m a.s.l.)L38J, and Tunka (in the Tunka Valley in Buryatia, Siberia, 675 m a.s.l.)||39l. A 
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Figure 4: Cosmic-ray energy spectrum measured with IT73 and other experiments operating in the same energy 
range. The spectra from the Pierre Auger Observatoryf32 J and Telescope Array[33 f are also shown. 
Finally, the H4a model of the cosmic-ray flux is overlaid {see text for detail). 



model assuming three populations of cosmic rays BOl (SNR component, high energy galactic 
component, and extra-galactic component) and termed as H4a is also shown. 

A first attempt to measure the mass composition of cosmic rays was performed with one 
month of data (constant snow coverage) of IceCube in its 2008 IT40/IC40 configuration llOI . 
A neural network was trained with Monte Carlo simulations of 5 primaries (proton, helium, 
oxygen, silicon, and iron). Measurements of the electromagnetic component of the air showers 
at the surface (S125) and the muon component in the ice (K70) are used to "teach" the network 
how to find the best fit to the primary energy and mass. A measurement of the cosmic- 
ray energy spectrum (Fig. |3j and composition (Fig. |5j at energies between 1 PeV and 30 PeV 
was determined. The energy resolution from the neural network ranges from 18-20% in the 
threshold region of this analysis (1.5 PeV) to 6-8% at 30 PeV, for an average resolution better 
than 14% over the full range of energies. The energy spectrum derived with this analysis is 
consistent with other IceTop results within systematic uncertainties. The mean logarithmic 
mass < InA > as a function of the primary energy indicates a strong increase in mass through 
the knee although the systematic uncertainties can greatly affect the measured composition in 
terms of absolute value of < InA >. 

IV. Cosmic-ray Anisotropy Measurements 

Although detection of 7 rays and neutrinos from individual galactic or extra-galactic sources 
of cosmic rays remains a method to probe the origin of cosmic rays, studies of anisotropy 
of cosmic-ray arrival directions are important to investigate the characteristics of cosmic-ray 
propagation in the local interstellar medium. 

In IceCube the anisotropy of cosmic -ray arrival directions can be measured in two different 
ways: using TeV muon events collected with the deep detector or using cosmic-ray air shower 
events triggering the surface array. The in-ice detector has a lower energy threshold than 
IceTop, which allows to investigate the anisotropy of cosmic rays at lower primary energies 
(down to 20 TeV) and larger zenith angles (up to 90°). This makes it possible for the deep 
detector to reach a higher sensitivity (about 6.3-10^'^ events/yr, anisotropy level S > 10~^) and 
scan small scale structures. On the other hand, IceTop presents a better energy resolution (20% 
at >300 TeV), although binning is limited by statistics, and potential for including composition 
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Figure 5: The mean logarithmic mass vs primary energy adapted from Ref. /TOf. The IT40/IC40 results are shown 
with (red) stars along with their statistical errors (solid red error bars) and systematic errors (shaded 
red region). Measurements from other experiments are also shown. Data points are compiled from 
Ref fMM- 



sensitivity. 

A dipole-like large scale anisotropy (amplitude of about 10^^) from few tens of TeV to 
about 100 TeV has been observed with in-ice data[42l |43l. This anisotropy is inconsistent, 
both in amplitude and phase, with the Comp ton-Getting prediction|44J, i.e. the apparent 
anisotropy caused by the relative motion between the Earth and sources of cosmic rays. A 
small scale anisotropy with significant structure at angular sizes between 10° and 30° has been 
also observed. This might uncover non-diffusive propagation effects or SNR connection and 
be a natural consequence of the stochastic nature of cosmic-ray galactic sources, in particular 
nearby and recent SNRs (<0.1-1 kpc)ll45ll46ll47l. 

The analysis of IceTop data confirms and complements the measurement of the large scale 
anisotropy also at PeV energies|48|. IceTop has an angular resolution of about 3° for energies 
above 100 TeV and zenith angles up to 40° when pure event geometry reconstruction is per- 
formed. The angular resolution degrades to 10° and above for zenith angles greater than 60°. 
Only events with zenith angles less than 60° are selected (1.4-10^ events/yr, S > 10^^). Monte 
Carlo studies indicate that the median primary cosmic-ray energy of IceTop data is 640 TeV, 
with 68% of the events between 200 TeV and 2,400 TeV. Events are classified as low (68% of 
events in 100-700 TeV, median energy of 400 TeV) and high energy events (68% of events in 
0.8-3.8 PeV, median energy of 2 PeV). 

Analysis of IceTop data reveals a deficit that confirms what was observed with the in-ice 
muon analysis at 400 TeV (compare lower left to upper right plot of Fig. [SJl. The amplitude 
of this deficit is about 210^'' and therefore larger than 10^'' observed with in-ice muons. 
However, these values are in agreement if the uncertainties are considered. Furthermore, above 
400 TeV there is indication of an increase in strength of the anisotropy. 

V. PeV 7 SEARCH 

High energy 7 rays (>1 TeV) have been observed from different galactic sources (SNRs, pul- 
sar wind nebulae, binary systems, the Galactic Center) and extra-galactic sources (starburst 
galaxies, active galactic nuclei, and objects containing supermassive black holes). At higher 
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Figure 6: Statistical significance sky maps obtained with IC79 in-ice (left panels) and IT73 IceTop data (right 
panels). The in-ice data sets have median energies of 20 TeV (upper) and 400 TeV (lower). The IceTop 
datasets have a median energy of 400 TeV (upper) and 2 PeV (lower). The angular binning or smoothing 
angle is 20°. 



energies (^100 TeV), extra-galactic photons are likely to interact with the cosmic microwave 
background radiation and radiation from infrared starlight from galaxies, producing — 
pairs. It is unknown whether galactic sources can emit 7 of energy >100 TeV. However, a 
guaranteed diffuse flux of 7 rays from interaction of cosmic rays with the interstellar medium 
and dense molecular clouds is also expected. 

IceCube can detect high energy 7 rays by vetoing on the in-ice component, i.e. searching 
for showers detected by IceTop where no in-ice activity is observed (muon poor showers) [i49i. 
A PeV 7-ray shower produces about 0.1 muons above 800 GeV The uncertainties due to the 
surface detector response and the muon rate production for photon showers are added in 
quadrature and return an overall systematic uncertainty in sensitivity of 18%. No detectable 
7-ray flux has been found by IceCube above 1.2 PeV with one year of data. The fraction of 
7 in cosmic rays is estimated to be less than 1.2-10^^ (90% cl) in the range 1.2-6.0 PeV and 
within 10° of the galactic plane. IceCube is also sensitive to localized sources, where galactic 
accelerators or dense targets for extra-galactic cosmic rays might be discovered. It is estimated 
that at about 1 PeV, IceCube can reveal fluxes as low as about 10^^^ cm^^s^^TeV^^ for point 
sources. 



VI. Laterally separated muons 

High energy muons (>1 TeV) are produced early in air showers and probe the initial shower 
development. One can distinguish between "conventional" muons, which come from pion 
and kaon decays, and "prompt" muons, which come from the decays of particles containing 
heavy quarks, mostly charm. The former are expected to dominate at TeV energies, the latter at 
higher but uncertain energies. Extending the earlier measurements performed by the MACRO 
experiment|51|, IceCube is capable to resolve muons of cosmic-ray primary interactions and 
high-energy secondary interactions that are laterally separated by 135 m up to 400 m from the 
shower core (Fig. EJ/ where the muon bundle is detected i52i|. The separation is due to high 
transverse momentum of 2-15 GeV/ c (corresponding to separations of 135-400 m) transferred 
to the muon by its parent. Above 2 GeV/ c, interactions can be described with perturbative 
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Figure 7: (a) Example of laterally separated muons (upper track) from the muon bundle (lower track) of a cosmic- 
ray shower reaching the in-ice detector. The separation is 400 m, the angular separation is 3.5°, and the 
time difference is 20 ns. (b) The separation between the laterally separated (LS) muon and the bundle 
track after applying all selection criteria for data, simulation with the Sibyll, DPMJET, and QGSJET 
interaction models, as well as double showers estimated from off-time data. 



quantum chromodynamics. However, the cosmic-ray hadronic interaction models used to 
simulate IceCube events do not reproduce the rates and the zenith angle distributions observed 
in data. 

VII. Solar and Heliospheric Physics, and GRB Search 

Due to the high altitude and the nearly zero geomagnetic cutoff at the South Pole, secondary 
particle spectra measured with IceTop retain a significant amount of information on the spec- 
tra of the primary particles. IceTop has already demonstrated the novel and unique ability to 
derive the energy spectrum of solar particles in the multi-GeV regime. The first event detected 
and studied by IceTop] 53 1 was the solar flare associated to the ground-level event of December 
13, 2006. IceTop DOMs can span thresholds from 1 pe to 30 pe corresponding to counting 
rates of about 8 kHz to about 1 kHz. By taking the differential rates at multiple thresholds, 
the spectrum of detected events can be studied. Furthermore, IceTop presents 50 times better 
sensitivity than conventional detectors used for solar and heliospheric physics such as neutron 
monitors. In addition to solar flares, the typical phenomena studied with IceTop are complex 
interplanetary disturbances! 55 1 and "Forbush Decreases"! 56 1. The latter are associated with 
strong shocks following coronal mass ejections that deplete cosmic rays from the region tra- 
versed by the Earth. On May 12, 2012 a significant solar flare was detected by IceTop and is 
being currently analyzed to determine the energy spectrum (Fig.!8]l. 

An ongoing study shows that IceTop can also reveal GRBs through detection of an overall 
increase of the counting rates observed in coincidence with the signal recorded on board of 
dedicated satellites. A flux of 7 rays with energies greater than 10 GeV and zenith angles 
up to about 20° can enhance the counting rate above the steady cosmic-ray background if it 
is greater than 10^^ erg cm^^. IceTop can detect GRBs whose emission extends up to a few 
100 GeV and occurring in a sky region not monitored by any other experiment. 
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Figure 8: Time profile of the count rate measured by groups of IceTop DOMs with decreasing thresholds (from 
MPE to SPE3), neutron monitors, and the GOES satellite^B^. 

VIII. Conclusions 

The IceCube Observatory is currently taking data in the second year after its completion in 
2010. Analysis performed with data of the detector in earlier stages of its deployment has 
been reviewed in this paper. On the one hand, events seen in coincidence by both the sur- 
face component, IceTop, and in-ice detectors offer precious information to investigate the mass 
composition of cosmic rays. On the other hand, analysis of IceTop events allows for precise 
measurements of the cosmic-ray energy spectrum with large sensitivity. The better under- 
standing of the systematic uncertainties along with the exquisitely large event rate collected 
will allow IceCube to significantly contribute to shed light on galactic cosmic rays in the near 
future. 
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